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Abstract 

CN . We propose to investigate infrared properties of gluon and ghost propagators re- 

lated to the so-called Gribov-Zwanziger confinement scenario, originally formulated 

■ for Landau and Coulomb gauges, for other gauges as well. We present results of our 
investigation of SU{2) lattice gauge theory in the maximally Abelian gauge (MAG), 
focusing on the behavior of propagators in the off-diagonal (i.e. non- Abelian) sector. 

(— I ' We also comment on our preliminary results for general linear covariant gauges, in 

particular for Feynman gauge. 



cn ; 1 Introduction 

^ . Important features of quark and gluon confinement in QCD are believed to be closely related 
Q ! to the behavior of gluon and ghost propagators in the infrared limit. One must notice, 
however, that the study of infrared properties of these propagators must be performed by 
nonperturbative methods and at a fixed gauge. The Gribov-Zwanziger confinement scenario 
[H [2] — proposed for Landau and Coulomb gauges — provides predictions for gluon and 
ghost propagators in the infrared limit, which may be tested by lattice simulations and 
by nonperturbative analytic methods such as Dyson-Schwinger equations. In particular, 
a suppressed infrared gluon propagator D{jP') is predicted, with -D(O) = 0. The latter 
statement implies maximal violation of reflection positivity for the gluons, a result that may 
be viewed as an indication of gluon confinement. (Note that it suffices to have violation of 
reflection positivity, not necessarily maximal.) At the same time, the infinite-volume limit 
favors gauge configurations on the boundary region known as the first Gribov horizon, 
where the smallest nonzero eigenvalue Xmin of the Faddeev-Popov matrix M. goes to zero. 
As a consequence, the ghost propagator GijP') — which is obtained from M.~^ — should 
be infrared-enhanced, introducing long-range effects in the theory. These, in turn, would 
be responsible for the color- confinement mechanism. 

The nonperturbative study of infrared propagators may be carried out from first prin- 
ciples in lattice simulations, taking into account that the true infrared behavior is however 
obtained only at large enough lattice volumes. Considerable effort has been dedicated to 
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investigations of the above predictions for the Landau gauge, considering very large lattice 
sizes (see e.g. [3]). The status of these studies is discussed in Here we propose to test 
similar predictions for the propagators as applied to the lattice implementation of other 
gauges, to try to gain a unified understanding of the mechanism of confinement and its 
manifestations. We consider the maximally Abelian gauge (MAG) and the linear covariant 
gauges, in particular Feynman gauge. 

In the case of the linear covariant gauges, which include and generalize Landau gauge, 
some studies suggest that the Gribov-Zwanziger confinement mechanism may apply to the 
complete class of such gauges [HE]- (A recent study of Dyson-Schwinger equations for 
Feynman gauge has been presented in [7].) On the other hand, for MAG, the usual confine- 
ment scenario is based on the concepts of Abelian dominance and of dual superconductivity 
[8]. Nevertheless, one might argue that a modified Gribov-Zwanziger scenario would likely 
hold in MAG for the non- Abelian directions in gauge-configuration space. A study of the 
Yang-Mills Lagrangian restricted to the (MAG) Gribov region by addition of a horizon 
function with Gribov parameter 7 has recently been carried out for SU{2) gauge theory 
in [9]. As pointed out in that reference and also by other groups, the infrared behavior of 
propagators in MAG may be modified by the presence of ghost and gluon condensates of 
mass dimension two. An example of such objects is the ghost condensate v [101 [H], related 
to the breakdown of a global SL{2, R) symmetry. This quantity is expected to modify 
the symmetric and anti-symmetric components of the (off-diagonal) ghost propagator. In 
particular, a nonzero value for v corresponds to nonzero anti-symmetric components of the 
ghost propagator. In Section[2]we present results of our lattice studies of pure SU{2) theory 
in MAG. (The implementation of gauge fixing for MAG on the lattice is straightforward.) 
We consider gluon and ghost propagators, the ghost condensate v mentioned above and the 
smallest eigenvalue of the Fadeev-Popov matrix. Our preliminary results have also been 
presented in |T2] and [13] . We note that the bounds recently introduced for studying gluon 
and ghost propagators on large lattices in Landau gauge [TH [15] may be written also for 
other gauges. 

Contrary to the case of MAG, the technical aspect of fixing the linear covariant gauges on 
the lattice is still not a settled issue. We comment on our recent proposals for gauge-fixing 
methods for these gauges in Section [31 

2 Infrared propagators in MAG 

g. [16]) by minimizing the 

(1) 

defined as 
-e,)[l-2{Ul{x)f] 

debcUi{x)Ul{x)]} , (2) 



On the lattice, for the SU(2) case, the MAG is obtained (see e.j 
functional 
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At any local minimum one has that the Faddeev-Popov matrix, 
Y.M^''{x,yh''{y) = J2r{x)[V,{x) + V,{x-e,)] + 2{r{x 



by 



2 Y.^\x - e,)[e,,U'^{x)Ul{x) + ^e. 



cd 



2 



is positive-definite. Here the color indices take values 1,2 and we follow the notation 
U,ix) = + ta'^U^^ix) and V,ix) = {U'^ix)f + {Uf,{x)f - {Ul{x)f - {Ul{x))\ 

where a"" are the 3 Pauli matrices. Notice that (as in Landau gauge [T7|) this matrix is 
symmetric under the simultaneous exchange of color and space-time indices. Using the 
relation U^{x) = exp [—iagoA^i^x)] one finds (in the formal continuum limit a — > 0) the 
standard continuum results [18] for the stationary conditions above and for M"*(x,y). 

We have considered four values of (3 (2.2, 2.3, 2.4, 2.512) and lattice volumes up to 40^. 
Data for a larger lattice, of volume 56^, have also been recently produced, but are not fully 
analyzed yet. We include data for the ghost propagator at this volume in Fig. [2] below, for 
comparison. 

Our results for the gluon propagators are in agreement with the study by Bornyakov et 
al. [16] : we see a clear suppression of the off-diagonal propagators compared to the diagonal 
(transverse) one, supporting Abelian dominance. We have fitted our data for the various 
gluon propagators (at all values of V up to 40"^ and for P = 2.2), obtaining the following 
behaviors. For D{p^) (transverse) diagonal, our data favor a Stingl-Gribov form 

^yp = ,u 2 I — I ' (3) 

with a mass m = \Ja/b ~ 0.72 GeV . Note that the above equation corresponds to 
a pair of complex-conjugate poles z and z*. We can thus write z = x + iy with x = 



b/{2c) ^ 0.32GeV^ and y = ^Ja/c-x"^ ^ OATGeV"^. Let us recall that in the case of 
a Gribov-like propagator these two poles are purely imaginary. For -D(p^) transverse off- 
diagonal our best fit is of Yukawa type, i.e. -D(p^) = l/(a + bp'^) , with a mass m = 

\J a/b ~ 0.97 GeV^ . Finally, the longitudinal off-diagonal gluon propagator is best fitted by 

= l/(a + bp^ + cp^) (i.e. also of Yukawa type) with a mass m = \J a/b ~ 1.25 GeV . 
As expected from Abelian dominance, the mass is larger in the off-diagonal case. 

In Fig. [1] (left) we show our data for the ghost propagator G'(p^), as a function of an 
improved momentum p (see Ref. [IS])- The data show little volume dependence at small 
p. (Note that, contrary to Landau gauge, here we can evaluate the ghost propagator at 
zero momentum.) We see no sign of an enhanced IR propagator. We have fitted our 
data (at (3 = 2.2), obtaining a behavior of the type ([3]) above with a = 0.45(1) GeV^^ , 
b = 1.1(3) , c = 0.73(30) GeV-^ , d = 2.1(9) GeV-^ . Thus, we see a Stingl-Gribov fit with 
mass m ^ 0.6 GeV and complex poles given by x ~ 0.75, y ~ 0.22. 

We next consider (see Fig. [H right) the smallest eigenvalue A^m of the Faddeev- Popov 
matrix. We have looked at Xmin for several lattice volumes and values of /5 as a function 
of 1/L. The data are fitted to a (l/L)'' with b = 1.6(1), showing that Xmin vanishes more 
slowly than {1/Ly (Laplacian). This may explain why we do not see a diverging ghost 
propagator at zero momentum even at rather large lattice volumes |15] . 

Following the analysis done in Landau gauge [17], we consider the anti-symmetric off- 
diagonal ghost propagator ( | eahG'^^ijP') /2 \ ) rescaled by L^/ cos (vrp^ a/L), as a function 
of the (unimproved) momentum p for all lattice volumes and (3 values considered. The 
data show nice scaling for all cases considered. The data at = 40^ and (3 = 2.2 can 
be fitted by $(p) = (a + bp/L'^){p^ + v^) with a = 0.002Q{7) GeV\ b = 32.6(7) Gel/-^ 
and = 1.7(1) GeV^ . We thus have a rather large ghost condensate v ~ 1.3 GeV^^ , but 
we cannot be sure that it survives in the infinite- volume limit, since the overall constant a 
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Figure 1: Left: plot of G'(p^) as a function of improved p for lattice volumes V = 16^, 24^, 
40'^ and /? = 2.2. Right: plot of the smallest eigenvalue of the Faddeev-Popov operator, as 
a function of the inverse linear size of the system. 



might be null. We can also fit data at several V^s and /3's for as a function of p and 

L (see Fig. ^ right). We obtain = (a + hp/L'^){p'^ + t;^) with a = 0.0033(6) GeV^ , 
b = 35.8{5) GeV~^ and = 1.87(8) GeV^^ . We note that the fit parameters change little 
with the (physical) lattice volume. In fact, data obtained recently for a larger lattice volume, 
56^, are seen to fall nicely on top of the fit done for the smaller volumes, as seen in Fig. [2] 
(right). 

We have also investigated possible effects of Gribov copies on our results, by consider- 
ing the difference between our standard gauge fixing (using the stochastic overrelaxation 
algorithm [20]) and the so-called smearing method [21j. The effects are found to be of the 
order of the statistical error. 



3 Linear covariant gauges 

As mentioned in the Introduction, gauge-fixing to linear covariant gauges (other than Lan- 
dau gauge) on the lattice is a challenge. More precisely, the gauge condition is given by 

d.A^ix) = A%x) (4) 

for real- valued functions A°-{x). As opposed to the case of Landau gauge — for which 
A°'{x) = and the gauge is fixed by minimizing a simple functional of the gauge-transformed 
links — in the general case no such functional exists ^22j. The solution to this problem 
presented in [22], based on the consideration of a modified gauge-fixing condition for the 
minimizing functional, may be affected by spurious minima and it leads to an altered form 
of the Faddeev-Popov matrix. We propose to consider a class of gauges on the lattice that 
coincides with the perturbative definition of linear covariant gauges in the formal continuum 
limit. Our method is based on a three-step process. Instead of minimizing a functional 
of A"'{x) directly, we first fix the gauge to Landau gauge, i.e. the transformed gauge fields 
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Figure 2: Left: plot of the quantity = L^/ cos (vrp^ a/L) ( | eabG'^^ijP') /2 | ) as a func- 

tion of p for lattice volumes V = 8^, 16^ 24^ 40^ and /3 = 2.2. Right: plot of ^{p'^) as a 
function of p and L. Data for a larger volume, 56^, are also included here. 

satisfy = 0. Then we determine a transformation 0^(a;) such that 



satisfies Eq. (|4]). Finally, we repeat the procedure for several functions A"(x) with a Gaus- 
sian distribution of width y^. The case ^ = 1 corresponds to Feynman gauge. The resulting 
distribution of dA"^{x) is shown for ^ = 1 in Fig. [3], in comparison with the original (Gaus- 
sian) distribution taken for A"(x). We see that the expexted distribution is fairly well 
reproduced. 

Our preliminary results were presented in [23j. We are currently investigating an alter- 
native method for fixing these gauges. 
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Figure 3: Distribution of dA"^{x) (solid line) compared with a Gaussian of width (dashed 
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